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1.  Introduction 


This  report  describes  the  mean  vertical  profile  of  the  refractivity  turbulence 
structure  constant  and  its  variability  with  season,  time  of  day,  and 
background  weather  conditions  at  White  Sands  Missile  Range  (WSMR),  NM. 

is  an  important  variable  for  studies  of  atmospheric  turbulence  and  for  the 
propagation  of  electromagnetic  waves  through  the  atmosphere. 
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2.  Data 


Data  are  from  the  50-MHz  radar  located  at  WSMR.  Details  of  the  radar  design 
and  operation,  along  with  examples  of  the  data  and  a  discussion  of  data  quality 
control  procedures,  are  given  by  Nastrom  and  Eaton  (1993a,  b).  The  basic 
period  of  record  available  for  this  study  was  January  1991  through  April  1994. 
Generally,  the  radar  operated  continuously,  providing  vertical  profiles  with 
150-m  resolution  from  approximately  5  to  20  km  every  3  min.  The  radar  has 
three  beams  (north,  south,  and  vertical),  and  each  beam  is  sampled  for  1  min 
in  a  continuous  rotation.  The  data  consisted  of  complete  profiles  of  the  first 
three  moments  of  each  Doppler  spectrum  (mean  Doppler  velocity,  spectral 
width,  and  returned  power  calibrated  as  C^^).  All  analyses  were  made  using 
hourly  mean  values  of  the  data  that  survived  quality  control  checks,  except  the 
hourly  standard  deviations  of  each  variable,  and  the  autocorrelation  functions 
of  log  are  based  on  3-min  observations.  Data  were  collected  during 
12,25 1  h  with  204,675  profiles  recorded.  The  number  of  hours  of  data  for  each 
month  is  given  in  appendix  A,  along  with  the  profile  of  the  percent  of  total  data 
available. 

Data  continuity  was  checked  by  making  time-series  plots  of  the  hourly  mean 
horizontal  and  vertical  windspeeds  and  standard  deviations,  the  hourly  mean 
spectral  width  (used  as  reported,  with  no  correction  for  effects  such  as  shear 
broadening),  and  the  hourly  mean  For  width  and  plots  of  the  monthly 
means  were  made  and  then  hourly  deviations  from  the  monthly  means  were 
plotted  to  permit  easy  comparisons  among  levels  in  the  vertical.  Examples  of 
the  plots  are  given  in  appendix  B;  a  complete  file  of  plots  for  January  1991 
through  April  1994  is  available  for  approved  requests.  These  plots  show  that 
the  radar  operated  for  only  6  h  per  day  (during  the  night)  from  October  27, 
1991,  to  June  2,  1992;  therefore,  these  months  will  not  be  used  in  studies  of 
diurnal  variations.  Some  brief  periods  had  anonymously  low  values;  these 
periods  were  removed  by  adding  a  requirement  during  data  checking  that  each 
profile  accepted  for  analyses  have  at  least  25  of  the  top  50  levels  pass  all 
quality  control  checks. 
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January  1994  through  March  1994  had  anomalous  values  most  of  the  time. 
For  example,  figure  1  shows  frequency  distributions  of  hourly  means  of  log 
at  5.6,  12.1,  and  17.0  km,  and  for  the  median  log  from  12  to  18  km  using 
data  from  the  9  winter  (December,  January,  and  February)  months  through 
December  1993.  The  distributions  appear  Gaussian,  as  expected  (past  studies 
such  as  Hufiiagel  1974;  Chadwick  and  Moran  1980;  and  Nastrom,  Gage,  and 
Ecklund  1986  have  found  that  the  distribution  of  appears  to  be  lognormal). 
Figure  2  shows  the  frequency  distributions  for  the  same  data  and  includes 
January  and  February  1994.  A  bimodal  distribution  is  evident  and  is  due  to  the 
spurious  values  from  1994.  Furthermore,  a  loose  electrical  connection  was 
found  in  the  radar  near  this  time.  Thus,  from  January  to  March  1994  will 
not  be  used  in  these  climatological  studies. 


WSMR  Log  C*  ot  5.6  km  N=2487  ,  rj  . 

(9  months:  gi;J-F-D  92:J-r-D  93;J-F-~0)  ^  ^  ^  N-2487 


#  XouanbSij  § 
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Figure  1.  Frequency  distribution  of  log  at  5.6,  12.1,  and  17.0  km,  and  for  the  median  from  12  to  18  km 
during  9  winter  months.  N  is  the  total  number  of  hourly  means. 


5.6  km  N  =  J412  WSMR  Log  oi  12.1  km  N  =  3412 

l-F-D  g3  J-r-D  g4:J-F)  (11  rrtonlhs:  91;J-F-D  92;J-F-D  93:J-F-D  94:J-F) 
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Figure  2.  Frequency  distribution  of  iog  at  5.6,  12.1,  and  17.0  km,  and  for  the  median  from  12  to  18  km 
during  11  winter  months.  The  very  low  values  of  are  spurious. 


3.  Seasonal  and  Interannual  Variations 


3.1  Mean  Horizontal  Winds 

Figure  3  shows  the  mean  profiles  of  zonal  u  and  meridional  v  wind  as  measured 
by  the  WSMR  50-MHz  radar.  The  number  of  hourly  means  available  at  7  km 
is  given  by  N  on  each  chart.  Error  bars  extend  ±1  standard  error  of  the  mean 

(SEM)  from  the  mean,  where  SEM  =  2a l\fN  and  a  is  the  standard  deviation 

of  the  N  values.  The  Jetstream  is  strongest  during  the  winter  (December, 
January,  and  February),  but  the  nose  of  the  jet  is  flat.  The  flat  nose  may  be 
understood  when  turbulence  is  very  low  in  the  center  of  the  jet,  and  thus  radar 
echoes  are  weak  in  the  center  of  the  jet  and  may  have  unacceptably  low  signal- 
to-noise  values.  Figure  4a  compares  the  seasonal  means  of  u  and  v.  The 
profiles  of  u  have  a  relatively  similar  shape.  The  shear  of  the  mean  wind  is 
about  the  same  during  all  seasons  (except  perhaps  summer)  in  the  stratosphere 
when  discussing  and  its  dependence  on  vertical  wind  shear.  Profiles  of  the 
seasonal  means  of  the  absolute  values  of  the  vertical  shear  of  zonal  and 
meridional  winds,  and  of  the  total  shear,  over  300-m  layers  based  on  individual 
profiles,  are  compared  in  figure  4b.  The  differences  among  seasons  are 
relatively  small  for  even  these  shears.  In  summer,  the  WSMR  radar 
observations  reach  into  the  stratospheric  easterlies.  Meridional  winds  are  from 
the  south  during  all  seasons  except  fall. 

3.2  Mean  Vertical  Winds 

Figure  5  shows  the  mean  vertical  velocity  w  and  the  mean  of  the  standard 
deviations  of  vertical  velocity  over  1-h  periods  (aj.  The  mean  w  is  downward 
at  several  centimeters  to  the  negative  first  power  in  the  troposphere,  and  is 
slightly  upward  in  the  lower  stratosphere.  Nastrom  and  VanZandt  (1994)  have 
presented  a  theory  predicting  that  the  observed  w  will  be  downward  when  there 
is  upward  energy  propagation  by  gravity  waves  due  to  nonuniform  radar 
reflectivities.  They  predict  that  the  magnitude  of  w  is  proportional  to  the 
gravity  wave  amplitude.  An  indicator  of  gravity  wave  amplitude  is  a^.  Mean 
profiles  of  w  and  are  compared  in  figure  6. 
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Figure  3.  Seasonal  mean  profiles  of  zonal  and  meridional  winds  at  WSMR. 


WSMR  SEASONAL  MEAN  W  AND  SIGMA-W  D-J-F  WSMR  SEASONAL  MEAN  W  AND  SIGMA- 
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WSMR  VERTICAL  BEAM  (m  s"')  WSMR  VERTICAL  BEAM  (m  S"') 

Figure  5.  Seasonal  mean  profiles  of  vertical  velocity  and  its  standard  deviation  over  1-h  periods  at  WSMR. 
(Note  that  is  the  standard  deviation  over  1  h.) 
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SICMA-W  (m  i“')  MEAN-W  (m  ■"’) 

Figure  6.  Seasonal  mean  profiles  of  vertical  velocity  and  its  standard  deviation  over  1-h  periods  at  WSMR 
with  all  seasons  plotted  together. 


In  the  troposphere,  the  largest  w  occurs  during  fall  and  winter.  The  largest 
also  occurs  during  fall  and  winter,  although  the  correlation  between  w  and  <j^ 
is  not  perfect.  The  correlation  between  w  and  is  not  perfect  because  trapped 
waves  add  to  but  make  no  contribution  to  w;  therefore,  the  relative 
proportion  of  trapped  and  propagating  gravity  waves  may  change  with  season. 

3.3  Spectral  Width 

Figure  7  shows  profiles  of  the  seasonal  mean  spectral  width  along  the  east  and 
north  oblique  beams.  Spectral  width  is  proportional  to  the  intensity  of 
small-scale  turbulence  and  the  eddy  dissipation  rate;  however,  as  discussed  by 
Hocking  (1985)  and  Gage  (1990),  other  factors  such  as  beam  broadening  and 
vertical  wind  shear  may  contribute  to  the  spectral  width  and  must  be  accounted 
for  before  using  width  as  a  turbulence  variable.  The  spectral  width  from  the 
vertical  beam  is  not  used  at  very  high  frequency  because  echoes  from  both 
turbulence  and  specular  reflections  are  expected.  In  figure  7,  note  that  the 
width  of  the  meridional  beam  is  greater  during  all  seasons  except  summer,  and 
that  the  difference  between  beams  is  usually  greatest  in  the  troposphere. 
Figure  8  compares  the  seasonal  means  from  the  zonal  beam;  the  greatest  widths 
are  seen  during  winter  and  spring.  In  the  troposphere  below  approximately 
7  km,  there  is  little  variation  with  season. 

3.4  Cn' 

Figure  9  shows  profiles  of  the  seasonal  mean  of  log  (because  the  logarithm 

of  is  normally  distributed,  its  use  permits  routine  statistical  procedures)  and 
compares  the  zonal  and  meridional  beams.  As  with  width,  the  value  of  the 
meridional  beam  is  greater  during  all  seasons;  however,  for  the  difference 
is  greatest  in  the  stratosphere  and  seems  to  be  approximately  the  same 
magnitude  during  every  season.  Figure  10  (right-hand  panel)  compares  the 
aimual  mean  profiles  for  the  two  beams;  the  zonal  beam’s  values  are  smaller: 
about  1  dB  in  the  troposphere  and  2  dB  in  the  stratosphere.  The  vertical  shapes 
of  the  mean  profiles  resemble  the  shape  of  the  vertical  profile  of  the  potential 
refractivity,  A^.  Af  is  proportional  to  the  humidity,  the  vertical  gradient  of 
humidity,  and  the  vertical  gradient  of  temperature  as  well  as  the  density  of  air. 
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Figure  7.  Seasonal  mean  profiles  of  spectral  width  at  WSMR. 


WSMR  SEASONAL  MEANS  OF  WIDTH 


Seasonal  mean  profiles  of  spectral  width  at  WSMR  with  all  seasons  plotted  together. 


WSMR  SEASONAL  MEAN  LOG  C 
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Figure  9.  Seasonal  mean  profiles  of  log  at  WSMR. 


WSMR  ANNUAL  MEAN  C§  (1/1991-4/1994) 


Annual  mean  of  log  (right  side)  and  seasonal  means  of  log  (left  side), 


(1) 


Local  values  of  are  given  by 


where 

a  =  a  constant  close  to  unity 

a'  =  the  ratio  of  the  eddy  diffusivities  and  is  near  unity 
Lq  =  the  outer  scale  of  turbulenee. 

hf  is  given  by 

=  (77*10"^— W—fl  +  _  1^500^  ^2) 

re  [6zi  T  )  IT  hz] 

where  the  symbols  have  their  usual  meanings  (Gage  1990).  Mean  values  of 
as  seen  by  the  radar  depend  on  the  fraction  of  the  sample  volume  that  is 
actively  turbulent. 

Variations  with  season  of  are  compared  in  figure  10  (left-hand  panel). 
Largest  values  of  log  are  in  the  troposphere  during  summer  due  to 
inereased  humidity.  In  the  troposphere,  the  single  largest  eontribution  to  the 
variability  of  Kf  comes  from  humidity  gradients  (Gage  1990;  Tsuda  et  al. 
1988).  During  the  other  seasons  and  in  the  stratosphere,  there  is  surprisingly 
relatively  little  variation  of  log  The  stratospheric  static  stability  does  not 
change  much  with  season  and,  as  seen  in  figure  4,  the  mean  wind  shear  (leading 
to  loeal  turbulence)  is  also  fairly  constant  with  season. 

Profiles  of  the  standard  deviation  of  log  are  given  in  figure  1 1.  The  curves 
labeled  HOURLY  are  the  standard  deviation  of  the  V  hourly  mean  values  of  log 
available  during  a  season.  The  hourly  curves  depict  the  low-frequency 
variations  of  log  and  contain  the  varianee  due  to  processes  such  as 
year-to-year  changes  and  diurnal  ehanges.  The  curves  labeled  HI  FREQ  are  the 
averages  of  the  standard  deviations  of  log  within  each  hour.  There  are 
usually  20  individual  observations  taken  in  an  hour.  Note  that  there  is  very 
little  change  among  seasons  of  the  high-frequency  variations. 
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Figure  11.  Seasonal  mean  profiles  of  standard  deviation  of  hourly  means  of  log  (HOURLY)  and  of 
standard  deviation  of  log  over  1  h  (HI  FREQ)  at  WSMR. 


Log  is  relatively  constant  with  height  in  the  stratosphere  in  figures  9  and 
10  and  knowing  the  median  value  (for  a  normal  distribution,  the  mean  and  the 
median  are  equivalent;  using  median  values  ensures  that  results  are  not  distorted 
by  an  unusual  outlier)  between  12  and  18  km  would  provide  a  great  deal  of 
information  on  the  profile.  Therefore,  the  median  value  between  12  and  18  km 
was  included  in  figure  1  in  the  fi-equency  distributions  during  winter. 
Frequency  distributions  of  log  at  5.1,  12.1,  and  17.0  km,  and  for  the 
median  from  12  to  18  km  during  the  other  seasons  are  given  in  figures  12 
through  14.  During  all  seasons,  the  distributions  are  broadest  at  5.6  km. 

Interannual  variations  of  the  seasonal  means  of  log  are  given  in  figure  15. 
The  interseasonal  variations  seen  in  figure  10  and  again  in  figure  15  are  clearly 
much  more  important  than  the  interannual  variations  seen  for  any  given  season 
in  figure  15,  because  every  year  the  mean  summer  profiles  in  the  troposphere 
are  the  largest  and  values  for  spring  and  winter  are  the  smallest.  Variations  are 
relatively  small  in  the  stratosphere.  There  is  no  apparent  trend  of  with  time 
at  any  height  at  WSMR,  in  contrast  to  the  year-to-year  trend  in  CO  reported  by 
Frisch  et  al.  (1990). 


WSMR  Log  Cj  at  5.6  km  N  =  2057  WSMR  Log  ot  12.1  km  N  =  2057 

(9  months:  91:M-A-M  92;M-M  95:M-A-M  94:A)  (9  months:  91;M-A-M  92:M-M  93:M-A-M  94:A) 
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18  km  during  9  spring  months. 


18  km  during  8  summer  months. 


WSMR  SEASONAL  MEANS  OE  LOG  CLi  1991-1993 


LOG  Gj^ 

Figure  15.  Comparison  of  individual  seasoual  means  to  show  interannual  variations  of  log 


4.  Autocorrelation  of  Log  Cjv,^ 

The  autocorrelation  functions  of  log  at  7,  12,  and  17  km  for  each  of  the 
three  radar  beams  are  given  in  figures  16, 17,  and  18  for  January,  February,  and 
July  1993,  respectively.  The  left-hand  panels  show  the  autocorrelation 
functions  fi-om  lag  0  to  about  lag  300  min  (5  h);  also  given  are  the  number  of 
data  pairs  at  lags  0,  3,  and  6  min.  The  right-hand  panel  shows  the  segment  of 
the  data  fi'om  the  left-hand  panel  for  lags  out  to  27  min  in  semilogarithmic 
coordinates.  These  results  suggest  the  autocorrelation  function  /-(J)  can  be 
modeled  fairly  well  as  the  sum  of  a  first-order  autoregressive  process  and  a 
random  process;  that  is, 

r(T)  = 

where  a<  \  and  v  <  0.  The  quantity  a  reflects  the  random  process.  The  decay 
rate  of  the  red-noise  process  is  v.  Taking  the  logarithm  of  equation  (3)  gives 

In  r(x)  =  In  a  +  VT  (^) 

which  shows  this  model  will  be  linear  in  semilogarithmic  coordinates.  Some 
of  the  curves  in  the  right-hand  panels  of  figures  16  through  18  appear  to 
decrease  rapidly  from  lag  3  min  to  lag  6  min  and  then  decrease  linearly  at 
higher  lags.  This  shape  implies  there  is  also  a  small  contribution  to  the 
autocorrelation  function  fi’om  a  moving-average  process,  quite  likely  due  to 
blobs  of  turbulence  that  overlap  firom  one  observation  to  the  next.  Surprisingly, 
there  is  little  apparent  difference  in  the  autocorrelation  fimctions  between  the 
troposphere  and  stratosphere.  Past  reports  of  persistent  layers  in  the 
stratosphere  as  reviewed  by  Gage  (1990)  might  have  led  one  to  expect  larger 
autocorrelation  values  in  the  stratosphere.  A  study  of  persistent  layers  at 
WSMR  is  currently  being  conducted. 

Using  equation  (4),  the  autocorrelation  function  can  be  integrated  to  estimate 
the  effective  time  T  between  independent  observations  using  the  relation 

T  =  2rr(z)dz.  (5) 

Jo 
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The  curves  in  figures  16  through  18  suggest  v  is  about  7.1  x  10'^  min"',  and  that 
0.5  <  a  <  0.9.  T  ranges  from  approximately  2  to  4  h  with  these  values.  This 
is  the  average  time  between  statistically  independent  observations  when 
employing  these  estimates.  The  time  before  another  observation  is  needed  may 
be  much  shorter  for  a  particular  situation  or  application. 


logplotc.g  Tue  Sep  27  09:22:59  1994  lagplotc.g  Tue  Sep  27  09:23:00  1994 

WSMR  reb93  LOG  AUTOCORRELATION  FUNCTION  WSMR  Feb93  LOG  C^  AUTOCORRELATION  FUNCTION 


Time  lag  (Minutes)  Time  lag  (Minutes) 

Figure  17.  Autocorrelation  function  of  individual  observations  of  log  during  February  1993  out  to  lag 
300  min  (left  side).  The  right  side  is  out  to  lag  27  min  and  in  semilogarithmic  coordinates. 


5.  Diurnal  Variations 


To  study  diurnal  changes  in  log  windspeed,  and  spectral  width,  the  data 
were  sorted  according  to  hour  of  the  day  during  each  season.  Time  series  of 
the  24-h  means  were  inspected  for  apparent  diurnal  cycles.  For  example, 
figures  19  through  22  show  the  diurnal  curves  for  beam  2  (the  meridional 
plane)  during  each  season.  In  each  chart,  seven  range  gates  have  been  averaged 
to  give  curves  for  layers  1.05  km  deep  for  speed,  width,  and  log  In  the 
lower  right  panel  of  each  chart,  20  range  gates  have  been  averaged  for  log 
to  give  layers  3  km  deep.  There  is  an  apparent  semidiurnal  variation  in  speed 
in  the  lower  stratosphere,  becoming  more  nearly  diurnal  at  the  upper  levels 
shown,  especially  during  the  fall  through  the  spring.  Width  and  log  have 
only  small  diurnal  changes  above  the  midtroposphere  during  all  seasons.  Width 
has  a  large  afternoon  maximum  in  the  troposphere,  especially  during  the 
summer. 

The  hourly  standard  deviation  of  vertical  velocity  has  a  relatively  large 
diurnal  variation.  Figures  23  through  26  show  the  diurnal  curves  of  over 
layers  1.05  and  3  km  deep  from  7.7  to  19.3  km  in  the  lower  panels.  Largest 
values  of  tend  to  be  found  in  the  afternoon,  near  22  UTC.  The  upper  panels 
show  w  and  width.  There  does  not  appear  to  be  a  large  change  of  the  diurnal 
cycle  in  with  altitude.  (The  spikes  of  seen  at  16  and  20  UTC  during 
spring  are  clearly  anomalous.)  No  pattern  of  diurnal  cycle  in  w  can  be  seen 
except  in  the  troposphere  during  spring.  Figure  27  shows  the  hourly  vertical 
profiles  of  for  each  season.  Clearly,  the  general  form  of  the  diurnal 
variation  in  does  not  change  much  with  height  in  a  given  season.  The 
diurnal  range  is  largest  during  summer.  Sato  (1990,  1992)  has  discussed  the 
diurnal  changes  in  vertical  wind  variability  at  the  MU  radar  in  Japan. 
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WINTER  Beam  =  2  SPEED  (1.05  km  loyer  averages)  WINTER  Beam=2  SPECTRAL  WIDTH  (1.05  km  loyer  averoges) 


Figure  19.  Diurnal  variations  of  windspeed,  spectral  width,  and  log  during  winter  on  the  beam  in  the 
meridional  plane. 


SPRING  Beam=2  SPEED  (1.05  km  layer  overages)  SPRING  Beam  =  2  SPECTRAL  WIDTH  (1.05  km  layer  overages) 


Figure  20.  Diurnal  variations  of  windspeed,  spectral  width,  and  log  during  spring  on  the  beam  in  the 
meridional  plane. 
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Figure  22.  Diurnal  variations  of  windspeed 
meridional  plane. 
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Hour  of  the  Day  Hour  of  the  Day 

Figure  23.  Diurnal  variations  of  windspeed,  spectral  width,  and  during  winter  on  the  beam  in  the 
meridional  plane. 


SPRING  Beam  =  3  SPEED  (1.05  km  layer  averages)  SPRING  Beam  =  3  SPECTRAL  WIDTH  (1.05  km  loyer  averoges) 


Hour  of  the  Day  Hour  of  the  Day 

Figure  24.  Diurnal  variations  of  windspeed,  spectral  width,  and  a,,  during  spring  oi 
meridional  plane. 
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6.  Changes  of  Log  with  Background  Weather  Conditions 


There  have  been  several  efforts  to  model  as  a  function  of  other,  more  easily 

measured  atmospheric  variables.  One  of  the  first  models  was  that  of  Hufiiagel 
(1974),  whose  independent  variable  was  mean  windspeed  aloft.  More 
sophisticated  models,  such  as  those  of  Wamock  and  VanZandt  (1985)  and 
Dewan  et  al.  (1993),  require  a  computer  to  produce  a  profile  of  and 
detailed  knowledge  of  the  vertical  profile  of  temperature,  humidity,  and  winds. 
The  goal  of  this  section  is  to  develop  simple  regression  models  that  could  be 
used  to  construct  an  estimated  profile  of  with  relatively  simple  and 
commonly  available  input  data  and  a  hand  calculator. 

Time-series  plots  of  Cn^  were  compared  with  other  variables  available  firom  the 
radar  (for  example,  horizontal  windspeed  and  direction,  and  standard  deviation 
of  the  vertical  velocity)  and  with  daily  weather  maps  at  the  surface  and  aloft. 
Also,  numerical  correlation  analyses  were  made  between  and  other  radar 
variables.  During  the  fall  through  the  spring,  the  tropospheric  responds 
strongly  with  the  airmass  changes  marked  by  the  passage  of  fi-onts.  Gossard 
(1977)  discusses  the  relationship  of  tropospheric  to  airmass  type.  A  good 

indicator  of  airmass  changes  is  the  meridional  windspeed  in  the  midtroposphere, 
near  500  hPa.  During  the  winter  season,  at  5.6  km  correlated  more  closely 
with  V  at  5.6  km  than  with  any  other  radar  variable,  including  spectral  width. 
This  relation  can  be  understood  as  the  winds  are  from  the  south  ahead  of  an 
eastward  moving  trough  at  500  hPa  and  fi-om  the  north  while  behind  it. 

The  close  correlation  of  log  and  v  is  illustrated  in  figure  28  using  data  firom 
January  1991.  The  upper  left  panel  of  figure  28  shows  time-series  plots  of 
meridional  windspeed  normalized  by  its  standard  deviation  over  the  233  data 
points  used  and  of  log  similarly  normalized.  The  correlation  is  0.407  and 
can  easily  be  detected  by  inspection.  The  upper  right  panel  shows  a  scatter  plot 
of  log  as  a  function  of  v,  along  with  the  linear  regression  line  for  this  data 
set.  The  regression  coefficients  for  all  data  during  winter  are  given  in  table  1 . 
The  lower  left  panel  gives  a  scatter  plot  of  log  C^^-U  versus  log  for 

January  1991.  That  the  two  variables  are  not  perfectly  correlated  may  result 
fi’om  different  sampling  volumes,  a  true  anisotropy  in  atmospheric  turbulence. 
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calibration  errors  in  one  or  both  radar  beam  systems,  or  other  factors  that  are 
beyond  the  scope  of  this  study.  The  lower  right  panel  of  figure  28  is  similar 
to  the  upper  right  panel,  except  for  log  as  a  function  of  meridional  wind. 

It  was  noted  earlier  that  Cj.,^  is  relatively  constant  through  the  lower 
stratosphere;  thus,  in  an  effort  to  keep  this  model  as  simple  as  possible,  the 
stratospheric  values  are  represented  by  the  median  value  fi-om  12  to  18  km. 
The  median  log  correlated  most  closely  with  the  standard  deviation  of 
vertical  velocity  in  the  stratosphere  at  12  to  18  km.  (Median  values  of 
are  summarized  in  figure  30  for  convenient  reference.)  The  physical  basis  for 
this  correlation  is  likely  that  increased  implies  increased  gravity  wave 
amplitudes,  which  in  turn  imply  more  wave  breaking  and  more  small-scale 
turbulence.  A  proxy  indicator  was  sought  because  is  not  a  commonly 
reported  variable.  The  median  was  found  to  correlate  closely  with  zonal 
windspeed  at  5.6  km  (near  500  hPa),  probably  because  the  mountain  ridges  near 
WSMR  run  generally  north-south,  so  the  zonal  wind  is  perpendicular  to  the 
ridges  and  thus  launches  vertically  propagating  gravity  waves  which  break  in 
the  stratosphere,  creating  turbulence  and  enhancing  Nastrom  and  Eaton 
(1993a)  also  found  strong  coupling  of  log  with  tropospheric  zonal 
windspeed.  Figure  29  shows  a  time-series  plot  of  zonal  windspeed  at  5.6  km 
and  the  median  log  at  12  to  18  km  for  a  period  in  spring  1991  (left  panel) 
and  a  scatter  plot  of  these  two  variables  (right  panel).  The  equation  for  the 
regression  line  in  the  right  panel  is  given  at  the  bottom  of  the  panel. 

Regression  coefficients  for  all  seasons  are  included  in  the  table  along  with  the 
correlation  coefficient  of  log  Cf.,^  with  the  independent  variable.  During  winter, 
the  correlation  of  log  (12  to  18  km)  with  zonal  windspeed  is  not  high.  If 
information  on  is  available,  it  should  be  used  instead  of  u,  since  it  has  a 
much  higher  correlation  (0.38).  The  proper  regression  equation  in  this  case  is 
~  -18.06  +  0.51*cj^y  (12  km).  The  value  of  at  5.6  km  does  not 
correlate  well  with  any  available  meteorological  variable  during  the  summer. 
Clearly,  in  the  troposphere  is  most  closely  related  to  the  level  of  humidity, 
and  it  is  likely  that  during  summer  humidity  is  generally  a  local  variable  not 
controlled  by  the  large-scale  flow. 


The  regression  model  implied  by  figures  28  and  29  and  the  above  discussion 
can  be  used  as  follows: 

•  Use  the  current  or  forecast  500  hPa  meridional  windspeed  to  estimate  log 
at  5.6  km. 

•  Use  ilie  current  or  forecast  500  hPa  zonal  windspeed  to  estimate  log  at  12 
to  18  km. 

•  Plot  the  estimates  on  a  graph  and  draw  a  line  from  the  value  at  12  km 
through  that  at  5.6  km. 

•  Extrapolate  the  estimates  below  5.6  km  to  the  top  of  the  planetary  boundary 
layer  (PBL).  This  model  is  for  the  free  troposphere  and  thus  not  appropriate 
for  the  PBL. 

Keep  in  mind  that  this  model  was  developed  using  data  from  WSMR  and  is 
likely  appropriate  only  for  similar  terrain  and  climatic  conditions.  Wamock 
et  al,  (1988)  and  Otten  and  Rose  (1985)  noted  differences  in  as  a  function 
of  underlying  terrain. 
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WSMR  Cj-U  ond  Meridional  Wind  ot  5.6  km  Jon  1991  WSMR  tog  C^-U  v3  Meridional  Wind  Speed  (5.6  km.  Jan  91) 


Figure  28.  Correlation  and  regression  relations  of  log  and  v  at  5.6  km 


Table  1.  Regression  coefficients  for  log  as  a  function  of  u  and  v  at  WSMR  (see  text) 
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WSMR  C8-U  (12-18km)  ond  Zonol  Wind  ot  5.6  km  M-A-M  1991  WSMR  log  Cj  (I2-I8km>  vs  Zonol  Wind  Speed  (5.6km)  M- 


Figure  29.  Correlation  and  regression  relations  of  log  (the  median  value  from  12  to  18  km)  with 
5.6  km. 
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7.  Summary  and  Conclusions 


A  relatively  very  large  data  set  was  collected  with  the  50-MHz  profiler  at 
WSMR.  The  available  data  span  the  period  January  1991  through  April  1994 
and  cover  the  height  range  from  approximately  5  to  20  km,  with  a  150-m 
vertical  resolution  and  a  3-min  time  resolution.  The  data  set  is  suitable  for 
climatological  studies  of  atmospheric  variables  such  as  windspeed  and  direction, 
vertical  wind  variability,  and  indicators  of  small-scale  turbulence  such  as 
spectral  width  and  This  report  has  emphasized  the  seasonal,  interannual, 
diurnal,  and  weather-related  changes  in  and  wind. 

Above  the  tropopause,  there  are  only  small  changes  with  season  of  Also, 
above  the  troposphere,  there  are  only  small  diurnal  changes  of  Cn^,  usually  less 
than  a  few  decibels. 

The  autocorrelation  function  of  can  be  modeled  as  the  sum  of  a  random 
process  and  a  first-order  autoregressive  process,  with  a  very  small  contribution 
from  a  first-order  moving-average  process.  The  shape  of  the  autocorrelation 
function  suggests  that  the  average  time  between  statistically  independent 
observations  of  ranges  from  about  2  to  4  h  in  both  the  troposphere  and 
stratosphere.  During  periods  of  active  weather  systems,  the  average  time  is 
much  less. 

Diurnal  variations  of  log  are  relatively  small.  Spectral  width  shows  a 
maximum  in  the  troposphere  during  local  afternoon,  as  does  the  standard 
deviation  of  vertical  velocity.  Mean  windspeeds  show  a  large  diurnal  variation 
in  the  stratosphere  during  summer. 

The  vertical  profile  of  is  related  to  temperature  and  humidity  gradients  and 
small-scale  turbulence  intensity.  As  a  result  of  these  relations,  in  the 
troposphere  is  correlated  with  airmass  changes  indicated  by  the  tropospheric 
meridional  windspeed,  whereas  in  the  stratosphere  is  correlated  with  gravity 
wave  activity  indicated  by  the  tropospheric  zonal  windspeed.  These  correlations 
were  used  in  developing  a  simple  regression  model  that  would  be  suitable  for 
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estimating  the  vertical  profile  of  with  only  simple,  commonly  available 
input  variables  and  a  hand  calculator. 

This  model  might  be  improved  by  sorting  the  data  according  to  another 
independent  variable(s)  such  as  temperature  or  static  stability.  No  rawinsonde 
data  were  used  in  this  study.  Also,  the  data  could  be  sorted  by  distance  fi’om 
the  tropopause  rather  than  by  geometric  altitude;  the  tropopause  is  a  dynamic 
variable  that  can  be  estimated  fi’om  the  backscattered  radar  power  (Gage  1990). 
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Appendix  A 


WSMR  Data  Inventory 


WSMR  Data  Inventory  1991—1994 
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Appendix  B 


Examples  of  Monthly  Charts 


Examples  of  the  charts  prepared  for  each  month,  January  1991  to  April  1994. 
These  charts  were  used  to  check  data  continuity,  for  comparison  with  standard 
weather  maps,  and  as  an  analyses  tool.  The  time-series  plots  are  of  hourly 
mean  values. 


WSMR  HORIZONTAL  WINDS  JAN91 
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Figure  B-1.  Horizontal  windspeeds  plotted  as  vectors  (estimated  from  the  oblique  radial  velocities  assumine 
the  vertical  velocity  is  zero). 
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Figure  B-2.  Standard  deviation  of  u  and  v  associated  with  the  means  in  the  chart  sho>yA<in  figure  B-1, 
plotted  as  vectors. 
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WSMR  MONTHLY  MEAN  WIDTH  JAN91 


Figure  B-6.  Vertical  profiles  of  the  monthly  means  for  spectral  width. 
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Deviations  of  spectral  width  from  the  monthly  means  in  the  chart  shown  in  figure  B-6 
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